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ABSTRACT 

We investigate the ability of the Space Infrared Telescope Facility's Infrared 
Array Camera to detect distant (z ~ 3) galaxies and measure their photometric 
redshifts. Our analysis shows that changing the original long wavelength filter 
specifications provides significant improvements in performance in this and 
other areas. 



Subject headings: cosmology: early universe — galaxies: evolution — infrared: 
galaxies — instrumentation: photometers — techniques: photometric 



1. Introduction 



On 1998 March 25, NASA Administrator Daniel Goldin approved the initiation of 
final design and construction for the Space Infrared Telescope Facility (SIRTF; Fanson et 
al. 1998). With launch planned for December 2001, SIRTF will complete NASA's family 
of Great Observatories. SIRTF will provide a large increase in sensitivity over previous 
missions across its 3-200 /im operating range, and over 75% of the observing time during 
its 2.5 year- minimum lifetime will be awarded to general investigators. A call for Legacy 
Proposals (large projects of both immediate scientific interest and lasting archival value, 
and with no proprietary data period) is planned for July of 2000. 

One of the defining scientific programs for SIRTF is the study of galaxies to z > 3 by 
means of deep surveys at 3-10 /im. This limit was selected because it is apparently beyond 
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the peak in the space density of luminous quasars (Schmidt, Schneider, & Gunn 1995). Not 
only will SIRTF's excellent sensitivity in this wavelength region allow such galaxies to be 
detected with the Infrared Array Camera (IRAC; Fazio et al. 1998), but the H~ opacity 
minimum at 1.6 //m (John 1988) is expected to be a major tool in photometric redshift 
determination at 1 ^ z <^ 5 (Wright, Eisenhardt & Fazio 1994), since it is a ubiquitous 
feature of stellar atmospheres. 

UV-bright examples of such galaxies have already been detected by means of the "UV 
dropout" technique (Steidel & Hamilton 1992, 1993; Steidel et al. 1996) and they play an 
important role in the overall star formation history of the Universe (Madau et al. 1996). By 
detecting galaxies on the strength of their UV emission, however, UV dropout samples are 
necessarily biased in favor of those with both active star formation and relatively modest 
extinction. Such samples will not reveal if there is an underlying population of galaxies 
which have already assembled the bulk of their stellar mass. In the absence of ongoing star 
formation, even massive galaxies will be too faint in the rest-frame ultraviolet to be picked 
up by optical surveys. The stellar mass already present at an early epoch constrains the 
star formation rate to that point, a quantity which is still uncertain due to the possibility of 
significant dust extinction. Steidel et al. (1998) suggest that "the onset of substantial star 
formation activity in galaxies occurs at z > 4.5", and estimate the net extinction correction 
to the star formation rate derived from UV dropouts to be a factor of about 5. Recent 
far-infrared detections of the cosmic background (Hauser et al. 1998), and far-infrared and 
submillimeter detections of field galaxies with ISO and SCUBA (Kawara et al. 1998; Puget 
et al. 1999; Barger et al. 1998; Hughes et al. 1998) also hint at a heavily dust enshrouded 
starburst population. Hence an accurate picture of the star formation history of the 
universe can only be determined by making an accurate census of all galaxies, not just 
star-forming ones with low extinction. Since the luminosity in the rest-frame near-infrared 
correlates linearly with mass (Gavazzi, Peirini & Boselli 1996) and is relatively unaffected 
by dust obscuration, this is clearly the spectral region in which to make such a census. 

In this paper we examine changing the specifications of the IRAC filters in order to 
optimize them for the study of z ~ 3 galaxies, including photometric redshift determination. 
IRAC will have two detector arrays made from indium antimonide (InSb) which is sensitive 
to A ^ 5 /im photons, and two made from arsenic-doped silicon (Si:As) which detects 
longer- wavelength light. All four arrays will have a 256 x 256 format, with 1'.'2 pixels, and 
will operate simultaneously. The InSb arrays will have filters centered at 3.63 /zm with 
20.4% bandwidth and at 4.53 /im with 23.3% bandwidth. We shall refer to these filters as L 
and M, respectively, because of their similarity to the filters used in ground-based facilities. 
The filters for the Si: As arrays have no ground-based analogs, since they cover wavelengths 
where the atmosphere is opaque, and so we refer to them as A and B, with A being the 
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bluer of the two. The L and A filters share one field and are separated by a dichroic filter, 
while the M and B filters similarly share a second nearly adjacent field. The original, or 
"baseline", specifications for A and B are 25% bandwidth filters centered at 6.3 and 8.0 /zm. 
Our analysis will involve comparing the performance of different pairs of filters, and we 
distinguish between these pairs by the use of subscripts, so for example the baseline filters 
described above are called A h and B h . We also introduce the notation [Ai : A 2 ] to describe 
a filter which extends from Ai /im to A2/xm. In this notation, the L filter is [3.26 : 4.00], 
M = [4.00 : 5.06], A h = [5.51 : 7.09], and B h = [7.00 : 9.00]. 



2. Assumptions 

The calculations in this paper all assume a total integration time of 10 hours per point 
(made up of many individual 200 second exposures), which allows a deep image of a large 
area of sky to be obtained in a sensible integration time. For such exposure times, IRAC 
requires background-limited operation and therefore we neglect read noise and dark current. 
Confusion noise will be discussed in §5.1 but is otherwise not considered here. As SIRTF 
is required to be diffraction limited at 6.5 /im, we adopt a detection area of twenty pixels, 
independent of wavelength, and assume that all the source flux is contained within this 
region. 



2.1. Telescope and instrument efficiency 

The fraction of photons incident upon SIRTF's 85 cm diameter aperture which reach 
the IRAC detectors is specified in the SIRTF Observatory Performance and Interface 
Control Document and the IRAC Instrument Performance Requirements Document^, and is 
about 50%, including filter transmission. The InSb quantum efficiency (QE) is assumed to 
be 80% over the entire wavelength range of the two InSb filters. The QE of the SiAs arrays 
is based on the average of laboratory measurements made on a number of flight-candidate 
arrays coated with anti-reflection coatings of either SiO or ZnS. Due to the high opacity of 
SiO beyond 8 /im, in the following analysis we assume that the ZnS coating will be used in 
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preference to the SiO coating if the filter bandpass extends beyond 8 /im. The net result 
of these assumptions is a QE of approximately 60-65% for the two Si:As filters. We also 
assume that the various filters considered here are opaque outside their nominal wavelength 
region. 



2.2. Background 



From the PICD, the minimum SIRTF background (which occurs at the ecliptic pole 
for wavelengths considered here), 6(A) is computed as 

6(A) = 5.5xl0" 14 J B A (5500)+4.2xl0" 8 J B A (278.5)+2.5xl0~ 4 (A/100/im) 2 J B A (17.5)+ J B A (2.73), 

where B\(T) is the Planck function for a temperature of TK. The first term is due to 
scattered sunlight from zodiacal dust, the second and third terms to thermal emission 
from zodiacal and Galactic dust, respectively, and the final term is the cosmic microwave 
background radiation. This function provides an excellent fit to data from the DIRBE 
instrument (Hauser 1996). 



2.3. Cosmology and galaxy models 



Our analysis is concerned with the detection of high redshift galaxies. We construct 
model spectra of such galaxies using the GISSEL96 spectral synthesis code (Bruzual & 
Chariot 1993, 1998; see also Chariot, Worthey & Bressan 1996). This code does not include 
polycyclic aromatic hydrocarbon (PAH) features at 3.3 /im and beyond, but at z = 3 these 
features are at observed wavelengths beyond the range considered here. We will return to 
the PAH features when we discuss confusion noise. 

We adopt a conservative (in the sense of making z ~ 3 galaxies faint) cosmology with 
H = 50kms _1 Mpc _1 and q = 0.1, corresponding to a present age of the Universe of 
16.5 Gyr (at z — 3 the Universe is 2.9 Gyr old). We assume a Salpeter (1955) IMF extending 
from 0.1-125 M , and normalize our galaxies to have M K = —24.63 (L*; Gardner et al. 
1997) at an age of 16.5 Gyr. We use four models to represent the limits of galaxy evolution. 
Model A is a maximally old model with solar metallicity formed in an instantaneous 
burst at z = oo and evolving passively thereafter. Model B is like A but with metallicity 
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only i solar. As we show in §3.3, these types of galaxies will not appear in UV dropout 
samples. Model C is a young solar metallicity galaxy formed with a constant star formation 
rate over the 100 Myr prior to whatever redshift is considered (i.e. a young nonevolving 
spectrum). Model D is like C, but again with Z = 0.02Z & . Models C and D are much more 
luminous than models A and B at z = 3, and are representative of the spectra expected 
for UV-dropout galaxies. The rest frame ultraviolet-optical-infrared spectra of these four 
models at z = 3 are shown in Figure |I[ 

The only feature common to all these spectra is a change in slope at 1.6 jum, due to 
the H~ opacity minimum, although this is barely discernable in Model D. It is in fact the 
only feature in the two low-metallicity models, which lack the prominent CO bandhead 
at 2.3 /^m and the effects of metal line blanketing at rest-frame optical wavelengths. This 
therefore makes it a suitable feature for use in photometric redshift determination, a 
technique whereby an approximate redshift can be inferred from observed colors as one or 
more spectral features pass through broad-band filters (e.g., Hogg et al. 1998). 



3. Sensitivity 

Since SIRTF's observations will be dominated by general observer programs, it is 
important to consider performance independent of any particular program, in addition to 
the high redshift galaxy observations which drive the filter recommendations made here. 
In fact, as we discuss in §6, the pair of filters which provide the optimum performance for 
A — B color selection and photometric redshift determination at z ~ 3 are not very suitable 
for other science programs. For this reason we will throughout the course of this paper be 
discussing the relative performance of three different Si:As filter pairs, which we designate 
"baseline" , "optimum" , and "preferred" and denote by subscripts b, o, and p, respectively. 
The baseline filters are the two 25% bandwidth filters centered at 6.3 and 8.0 fim, the 
optimum filters provide the most accurate A — B color selection and photometric redshift 
determination (see §4) at z ~ 3, and the preferred filters give slightly poorer performance 
in these areas, but are expected to be much more useful for other scientific programs. The 
optimum filters are A Q = [5.06 : 6.00] and B Q = [6.30 : 10.20], and the preferred filters are 
A p = [5.06 : 6.50] and B p = [6.50 : 9.50]. 



3.1. Power law spectrum source 
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Figure ^ is a contour plot of IRAC's background limited 5a sensitivity in 10 hours of 
integration as a function of filter cuton and cutoff wavelength for a source which has a 
flat spectrum in S v . Since the contour lines are fairly vertical, it is evident that the cuton 
wavelength is more important than the cutoff. It is also clear from this figure that the best 
sensitivity is achieved for a very blue filter whose short wavelength cuton will be constrained 
by the dichroic response. Within the parameter space plotted, a [4.2 : 5.7] filter is most 
sensitive, with a 5a 10 hr limiting flux of 1.1 /xJy However, this filter overlaps the InSb M 
filter almost entirely, and since InSb has a higher quantum efficiency there is little point in 
duplicating it. We therefore restrict our analysis to filters with A cut0 n > 5.06 /xm (the cutoff 
of the M filter). The second column in Table [I] lists the sensitivity achieved by various 
IRAC filters for a flat spectrum source. The most sensitive Si:As filter is [5.06 : 7.4] but, as 
remarked above, the sensitivity is largely independent of the cutoff wavelength, dropping 
only 1% if we reduce the cutoff to 7/xm or raise it to 8.3 /im. 

The optimum and preferred filters have a ~ 20% gain in sensitivity over the baseline 
specifications for a flat spectrum source. Figure ^| shows that the improvement holds for 
a fairly large range of spectral shapes, with the preferred A p filter only losing out to the 
baseline A h filter for spectra redder than a ^ 3 (S u oc v~ a \ due to the longer wavelength of 
the Ah filter. 



3.2. High redshift galaxies 

Figure |] plots the SNR achievable in 10 hours through Si:As filters for the four galaxy 
models at z — 3. Again the contours are fairly vertical showing the cuton wavelength is 
important and the SNR achievable is very insensitive to the cutoff. The best SNR at z = 3 
for a Si:As filter is again with the bluest wavelength cuton allowed (5.06 /im) and with a 
cutoff of 7.4, 6.1, 7.2, and 6. 7 /mi for models A, B, C, and D, respectively. Table p] lists the 
signal-to-noise ratios achievable at z = 3 for the various Si:As filter choices. 

The general conclusions thus far are that substantial improvements in sensitivity can 
be achieved over the baseline filter specifications by using a bluer A filter and a broader B 
filter. We show that this applies over a broad range of redshift in Figures [5] and || which 
plot the signal-to-noise ratio achievable in 10 hours for the three different filter pairs. For 
completeness, we show the SNR achievable through the two InSb filters as a function of 
redshift in Figure |7|. The higher quantum efficiency of InSb and the lower background at 
short wavelengths produce an order of magnitude increase over what is achievable with the 
two Si:As arrays, but the need for the Si:As arrays will become clear in §4. 
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3.3. The SIRTF advantage 



The advantage in using SIRTF to obtain an effectively mass-limited sample of 
high-redshift galaxies has already been stated. Here, we compare the sensitivities calculated 
above to the limits obtained, or obtainable, by other methods. 

In Table |2], we first list the limits obtained, or obtainable, by deep exposures in the 
optical-infrared region. In Figure [8], we present the observed flux of a solar metallicity z = 3 
L* galaxy, which was formed in an instantaneous burst, through three different HST filters 
and our preferred IRAC filters, as a function of age. These fluxes correspond to model A 
where they intersect the right axis. From this, it is clear that optical selection methods, such 
as the "UV dropout" technique, can only find galaxies which have undergone substantial 
star formation since z = 4 (i.e., similar to models C and D), as passive evolution since 
then will cause their UV flux to drop below a plausible detection threshold. Near-infrared 
observations fare rather better, although it should be noted that the NICMOS HDF 
catalogue suffers significant incompleteness even at high signal-to-noise ratios (the 80% 
completeness limit is tuab ~ 26.7 in both filters; Thompson et al. 1998). This F160W limit 
is sensitive enough to detect old, sub-L* galaxies, and is comparable to the IRAC filters in 
this respect. However, the F110W filter, by virtue of its sampling a region below the 4000 A 
break, is far less sensitive, and obviously there is limited science which can be done from 
a single detection in one filter. In addition, the short rest-frame wavelengths which these 
filters probe are strongly affected by dust obscuration. IRAC will therefore be the first 
instrument capable of detecting distant, evolved L* (and possibly even sub-L*) galaxies, 
if they exist, at more than one wavelength. As a comparison, we note that the deepest 
observations obtained by ISOCAM with the very wide LW-2 filter (5.0-8.5 /zm) produced 
a much higher 5a limit of approximately 30 /iJy in 6 hours of integration (13 hours of 
observation; Taniguchi et al. 1997). 



4. Photometric redshift determination at z ~ 3 



In this section we present the analysis which defines the "optimum" filters and leads to 
our specifications for the "preferred" filters. 

Two factors are important in deciding upon the best filters for photometric redshift 
determination. First, it is necessary to achieve a high signal-to-noise ratio in both filters, 
to allow the color to be determined accurately. This was discussed in the previous section. 
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Secondly, this color must be a strong function of redshift, since the uncertainty in redshift, 
Az, is related to the uncertainty in the color, Ar (the color defined in the sense T = A — B) 
by 

Az = (dT/dz)- 1 AT. (1) 

Figure || plots the L — B, M — B, and A — B colors as a function of redshift for the four 
models. Magnitudes in each of the filters have been determined with respect to Vega (which 
we model as a 9400 K blackbody with a 2.2 /im flux of 657 Jy). It should be noted that at 
z <; 2, the 3.3 /im PAH feature may affect the results since it will lie in the bandpass of filter 
B. It is apparent that the L — M color is sensitive to redshift in the range 1 < z < 2 (see the 
divergence between the L — B and M — B loci over this redshift range), due to the 1.6 /im 
bump moving through the filters. At higher redshifts the Si:As filters become increasingly 
important for measuring photometric redshifts. At z > 2, the increased separation between 
the preferred and optimum A and B filters, compared to the baseline filters, causes them 
to produce larger color variations. Coupled with the increased sensitivity that these filters 
provide, it is clear that they will be able to produce more accurate photometric redshifts. 
We now quantify the level of this improvement. 

The goal of photometric redshift determination is to minimize the uncertainty in the 
derived redshift, given by Equation [I]. We make a few approximations to ease the calculation 
of this function. First, we write dT/dz ~ T(3.5) — T(2.5), and make use of the fact that for 
well-detected (^ 5a) sources, the uncertainty in the measured magnitude Am « (SNR) -1 , 
where SNR is the signal-to-noise ratio. We approximate the SNR (which varies slightly 
with redshift) as the geometric mean of the actual SNR at z = 2.5 and z = 3.5. In light of 
these rather crude approximations, we do not claim to be deriving accurate values of Az, 
but instead present its reciprocal, which we call / (by plotting the reciprocal, we also avoid 
the infinities which arise when the color change between z = 2.5 and z = 3.5 is zero). For 
any two filters, X and Y, f is therefore defined as 

f(X Y) - r(3.5)-T(2.5) (2 ) 

\c x (3.5)c x (2.5) ^ c Y (3.5)c Y (2.5)J 

where 6j is the total number of background counts and Cj(z) the total number of object 
counts in filter i for a galaxy at redshift z. Large values of / are desired, and for photometric 
redshifts accurate to 10% at z ~ 3 (Az <^ 0.3), / ^ 3 is necessary. 

We wish to find the specifications for the two Si:As filters that maximize f(A,B). 
Unfortunately, since each filter needs two parameters to specify it, it is impossible to 
graphically present the results of this four-parameter optimization in a simple form. We 
therefore first consider the value of / for a single Si:As filter and the InSb M filter (the 



- 9- 



results are broadly similar if the L filter is used instead of M). Figure 10 is a contour plot 
of / for the M filter and an arbitrary Si:As filter whose specifications are given by the two 
axes. This figure favors a Si:As filter centered at ~ 8 fim with ~ 50% bandwidth. Note 
that the contours are predominantly horizontal here, so that contrary to the SNR case, the 
long wavelength cutoff is the most important parameter. A broad, long wavelength filter is 
optimum for determining photometric redshifts at z ~ 3 even though the SNR it produces 
is not maximal. 

We find the solution to the four-parameter optimization by performing a grid search 
throughout the available parameter space. We find that, irrespective of the particular galaxy 
model, the A filter always has the bluest cuton allowed (5.06 /xm). This is to be expected, 
since a blue A filter produces both high signal-to-noise (Figures and |j) and a large color 
change over the redshift range of interest. The other parameters vary slightly between the 
four galaxy models, but an average solution is A = [5.06 : 6.0] and B = [6.3 : 10.2]. This 
is the how we define the "optimum" filter choice referred to throughout this paper. The 
optimum B filter has a similar specification regardless of whether the color is formed using 
L, M, or A as the second filter. 

Table [3] lists the values of / achieved for various pairs of filters. In terms of the value 
of this function, the optimum filters can produce an increase of as much as 60% in the 
accuracy of the photometric redshifts derived from the A and B filters, and a more modest 
(but still substantial) 20% improvement in accuracy when one of the InSb filters is used in 
conjunction with the B filter. 

As an illustration of how photometric redshifts might be derived from IRAC data, 
in Figure [XT| we present a L — M v s A p — B p color-color plot for the four galaxy models 
in the redshift range 1 < z < 5. It can be seen that the color in the two Si:As filters is 
generally able to provide an excellent measurement of the galaxy redshift for z ^ 2. For 
1 ^ z ^ 2, the L — M color provides most of the photometric redshift signal. Only Model D 
has a very limited range of colors which might hamper analysis, since we are observing the 
Rayleigh- Jeans tail of a recent starburst with only weak metal line blanketing; however, 
this is exactly the sort of UV-bright galaxy which would be detected in surveys for UV 
dropouts, and so this does not pose a problem. 



5. Other considerations 



Finally, we consider the impact of filter changes on two important aspects of distant 
galaxy detection which we have thus far ignored: confusion noise, which will ultimately 
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limit SIRTF's sensitivity, and dust, which is not included in the Bruzual & Chariot spectral 
synthesis models. 



5.1. Confusion noise 

Confusion noise is expected to limit IRAC sensitivity to ~ 0.5 /iJy (la; Franceschini 
et al. 1991), although this number is quite uncertain and model-dependent. The confusion 
noise also depends on the source density and the beam size, both of which are wavelength 
dependent. Of particular concern is how much the density of confusing sources will increase 
over the baseline filter with a new B filter, due to a longer cutoff wavelength. The longer 
cutoff may significantly increase confusion because it includes the strong, broad 7.7 /im 
PAH feature (e.g. Uchida, Sellgren & Werner 1998) out to z ~ 0.5, a large enough distance 
that substantial volume is sampled and evolving dwarf galaxies may dominate the counts. 

Franceschini et al. (1991) estimate the density of sources at 6.7/im, and their prediction 
of ~ 1.2arcmin~ 2 above a flux level of 30/iJy is consistent with Taniguchi et al.'s (1997) 
detection of 15 such sources in a 3' x 3' field. At microJansky flux levels, most sources 
will be at z ^ 1, so we use the spectrum S v oc z/~ - 3 displayed at rest wavelengths ~ 3 /iin 
to convert a limiting 8 /xm flux to a 6.7/im flux. Our best estimate of the total number 
of sources above 0.5 /iJy at 8/xm is therefore ~ lOOarcmin -2 , which is ~ 0.8 sources 
per 20 pixel detection area. We estimate the number of PAH-emitting galaxies from the 
cirrus/photodissociation region spectrum and local 15 /im luminosity function of Xu et al. 
(1998). Even with very strong ((1 + z) 4 ) evolution, the number of sources with z < 0.5 
above 0.5 /iJy is no more than lOarcmin -2 , a small fraction of the total confusing source 
density. 

We also note that the longer cutoff wavelength of the new B p filter does not produce 
a significant increase in the diffraction-limited beam size for a range of spectral slopes, 
compared to the B b filter (Gautier, private communication). The confusion noise will 
therefore not be increased and the overall performance of the filter will not be affected. 

Finally, we note that the confusion noise limit of 0.5 /iJy would limit the ultimate 
values obtainable for f(A p , B p ) (as defined in Equation 2) to 3.2, 1.5, 34, and 3.7 for Models 
A, B, C, and D, respectively. In the confusion-limited case, where the photometric accuracy 
is largely independent of the filter choice, it can be argued that one should attempt to 
maximize the product of the photometric redshift accuracy and the sky area which can be 
surveyed in a given time. This produces slightly different specifications for the optimum B 



- 11 - 



filter ([6.0:10.1] instead of [6.3:10.2]) although the arguments of §6 lead to the same final 
choice of 'preferred' filters. 



5.2. Dusty galaxies 



Nearby starburst galaxies, e.g. M 82 and Arp 220, are very dusty and have red 
near-infrared spectra, as well as strong PAH emission features. Our choice of a bluer A filter 
might therefore impair our ability to detect similar objects at high redshift, which would 
be cosmo logically very interesting. Figure || shows that for spectra redder than a ^ 3, the 
A\y filter is more sensitive than A p . In Figure we show the sensitivities of the three 
different filter pairs to these two starburst galaxies as a function of redshift (the spectrum 
of M 82 has been scaled up by a factor of 10 in luminosity in this figure). The spectral 
energy distributions required to make these figures were produced with photometric data 
culled from a variety of sources, and aperture matching using curve-of-growth analysis was 
required. This is likely to produce errors where the spectral region sampled by a given filter 
includes the transition between different datasets (e.g. optical and near-infrared), and the 
signal-to-noise ratios (both absolute and relative) should probably not be trusted to more 
than ~ 10%. Given these caveats, there is little to choose between the different A filters, 
but the sensitivity improvement offered by the B p filter is unambiguous. 



6. Conclusions: the "preferred" filters 



Although the broad B Q filter is supported by sensitivity calculations and does not suffer 
significantly from increased confusion noise, a somewhat narrower filter makes sense for 
other reasons. The almost 50% bandwidth of the B Q filter will make photometry difficult 
and may compromise other uses of this filter. In particular, the strong silicate absorption 
feature at 9.7/im falls in the B Q bandpass. We therefore recommend a shorter cutoff 
wavelength, and propose B = [6.5 : 9.5] (a 37% bandpass) as our "preferred" B filter, B p . 

Conversely, the optimum A p = [5.06 : 6.0] choice is only a 17% bandwidth, and does 
not produce optimum SNR (see Table |l|). It does however avoid the PAH feature at 6.2 /xm, 
offering the possibility of using filters A and B as a PAH diagnostic. However the M filter 
also acts as a good continuum filter for this problem, with better SNR and the advantage 
that it shares the same field as the B filter. Hence we recommend a broader A filter filling 
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the wavelength regime between M and B. Our preferred A filter, A p , is then [5.06 : 6.5], 
with a 25% bandpass. 

These preferred filters reduce the value of / by about 10%, compared to the optimum 
filters, but we feel that this is outweighed by the gain in general scientific usefulness. In any 
case, they still offer a significant improvement over the baseline 6.3 and 8.0 filters. These 
specifications have now been adopted for the IRAC instrument. 
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Fig. 1. — The rest-frame spectra at z = 3 of the four galaxy models. The solid horizontal 
bars indicate the preferred filter set, while the dashed bars show the optimum Si:As filters 
and the dotted bars the baseline Si: As filters. 
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Filter cuton (/xm) 



Fig. 2. — Contours of limiting flux in fiJj (5<r, 10 hours) as a function of filter cuton and 
cutoff wavelength for a source with a flat spectrum (in S„). The contour range is from 
0.6 to 5.0yuJy in steps of 0.2 /iJy, with heavier lines at 1,2,3,4,5 //Jy. The discontinuity at 
^cutoff = 8/xm is due to the change in coating from SiO to ZnS. Squares mark the preferred 
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Spectral index a Spectral index a 

Fig. 3. — Limiting sensitivity (5er, 10 hours) for the preferred (solid line), optimum (dashed 
line), and baseline (dotted line) filters for power laws of the form S u oc v~ a . Left: Filter A. 
Right: Filter B. 
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Fig. 4. — Contours of signal-to-noise ratio as a function of filter cuton and cutoff in a 10 hour 
integration on the four galaxy models at z — 3. Symbols mark the filter specifications as in 
Figure ^|. Model A is contoured from 2 to 38 in steps of 2; Model B from 2 to 34 in steps of 
2; Model C from 20 to 480 in steps of 20; and Model D from 10 to 210 in steps of 10. Heavy 
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Fig. 5. — Signal-to-noise ratio (in 10 hours) for the three different A filters, as a function of 
redshift for the four galaxy models. Solid line: preferred (A p ); dashed line: optimum (A Q ); 
dotted line: baseline (A h ). 
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Fig. 6- 



As Figure ||, but for the three B filters. 
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Fig. 7. — As Figure ||, but for the InSb L (solid line) and M (dashed line) filters. 
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Age (Gyr) 



Fig. 8. — Observed flux for an L* galaxy at function of age (or formation redshift, 

z f ), through the filters ffi>T/WFPC2 F606W, HST /NICMOS F110W, #5T/NICMOS 
F160W, and ground-based K (all as labeled), and the four SIRTF /IRAC filters L (dotted 
line), M (dot-dash line), A p (dashed line), and B p (solid line). The sensitivities of the filters 
are listed in Table |2] and the one-sided error bars indicate the amount of dimming caused by 

Ay = 1. 
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Fig. 9. — Colours as a function of redshift for the four galaxy models. Solid lines: preferred 
filters; dashed lines: optimum filters; dotted lines: baseline filters. 
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Fig. 10. — Contours of / (as denned in the text) for the InSb M filter and an arbitrary (ZnS 
coated) Si:As filter as defined by the axes of the plot, for each of the four galaxy models. 
Model A is contoured from —0.4 to 6.2 in steps of 0.2; Model B from 0.0 to 3.1 in steps of 
0.1; Model C from -16 to 58 in steps of 2; and Model D from -1.2 to 6.6 in steps of 0.2. 
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Fig. 11. — Loci of the four galaxy models discussed in the text in the L — M vs A p — B p 
color-color diagram. The locations of the models at z = 1,2,3,4,5 are indicated. 
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Fig. 12. — Signal-to-noise ratio (in 10 hours) for redshifted versions of the starburst galaxies 
M 82 and Arp 220. M 82 has been scaled up by a factor of 10 in luminosity. Solid lines: 
preferred filters; dashed lines: optimum filters; dotted lines: baseline filters. 
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Table 1: Signal-to- noise ratios achievable with IRAC in 10 hours. 



SNR for Model 



Filter 


10/iJy 


A 


B 


C 


D 


Best filter a 


55.9 


26.9 


20.0 


309 


105 


^cutoff 


[6.9] 


[7.1] 


[6.7] 


[6.9] 


[6.4] 


L 


232.2 


70.9 


82.4 


1104 


649 


M 


136.7 


53.4 


49.5 


688 


321 


A h 


35.4 


21.8 


15.4 


247 


77 


A 


52.6 


23.6 


18.8 


281 


103 


A P 


55.5 


25.9 


19.9 


304 


105 


B h 


23.6 


10.3 


6.4 


112 


30 


B 


31.4 


12.8 


8.3 


139 


39 


B P 


28.9 


12.6 


7.9 


136 


37 



a The "best filter" is the most sensitive for a lO^Jy fiat spectrum source and the four galaxy models at z = 3, 
with the constraint that no filter extends bluer than 5.06 /im. All such filters cut on at 5.06 /xm and extend 
to the cutoff wavelength listed. 



-28- 



Table 2: Sensitivities (5a, 10 hours) of instrument /filter combinations. 



Instrument / filter 




Su (jMJy) 


HST WFPC2/F606W a 


28.4 


0.016 


HST NICMOS/F110W a 


27.2 


0.047 


HST NICMOS/F160W a 


27.5 


0.037 


Subaru IRCS/K b 


24.4 


0.63 


ISO CAM/LW-2 C 


20.4 


24.7 


SIRTF IRAC/L 


25.6 


0.21 


SI RTF IRAC/M 


25.0 


0.36 


SIRTF IRAC/A P 


24.0 


0.90 


SIRTF IRAC/^p 


23.3 


1.73 



a Scaled from observations of the Hubble Deep Field (Williams et al. 1996; Thompson ct al. 1998) 

fc A description of IRCS can be found in Tokunaga et al. (1998) 

c Scaled from Taniguchi et al.'s (1997) observations of the the Lockman Hole 
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Table 3: Values of / for various filter pairs. 





f for Model 


Filters ABC 


D 


Baseline filters 


L 


B h 6.4 3.1 57 


7.6 


M 


B h 5.4 2.6 52 


5.7 


A h 


B h 4.1 1.6 41 


4.3 


Optimum filters 


L 


B Q 7.6 3.7 66 


9.1 


M 


B Q 6.2 3.1 59 


6.6 


A 


B 6.0 2.8 63 


6.7 


Preferred filters 


L 


B p 7.3 3.5 65 


8.8 


M 


B p 6.0 3.0 58 


6.4 




B p 5.6 2.5 58 


6.1 


A 


B p 5.9 2.7 62 


6.6 



